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There exist strong evidence supporting the co-evolution of central supermassive black holes and 
their host galaxies. It is however still unclear what the exact role of nuclear activity, in the 
form of accretion onto these supermassive black holes, in this co-evolution is. We use a rich 
multi-wavelength dataset available for the North Ecliptic Pole field, most notably surveyed by the 
AKARI satellite infrared telescope to study the host galaxy properties of AGN. In particular we 
are interested in investigating star-formation in the host galaxies of radio-AGN and the putative 
radio feedback mechanism, potentially responsible for the eventual quenching of star-formation. 
Using both broadband SED modeling and optical spectroscopy, we simultaneously study the nu- 
clear and host galaxy components of our sources, as a function of their radio luminosity, bolo- 
metric luminosity, and radio-loudness. Here we present preliminary results concerning the AGN 
content of the radio sources in this field, while offering tentative evidence that jets are inefficient 
star-formation quenchers, except in their most powerful state. 
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1. Introduction 

The discovery of a number of scaling relations (e.g., Ferrarese & Merritt 2000, Magorrian 
et al. 1998, Kormendy & Richstone 1995) connecting nuclear properties of galaxies to their global 
characteristics has led to an ongoing debate about the possible physical processes that give rise to 
this connection. 

One such candidate physical process is feedback. The presence of activity within the nucleus of the 
galaxy (AGN) can give rise to outflows that may affect the whole galaxy. In particular signatures of 
star-formation quenching should be found in the host galaxies of radio-AGN, as radio-loud active 
galaxies are characterized by powerful, well-coUimated outflows. 

AGN radio-jets are known to be able to deposit large quantities of mechanical energy in their sur- 
roundings (e.g., McNamara et al. 2005). The role of radio-loud AGN and their jets with respect 
to star-formation has been studied intensively (e.g., Bicknell et al. 2000, Croton et al. 2006, Best 
et al. 2007; Best & Heckman 2012, Kalfountzou et al. 2012). Results have been ambiguous, with 
evidence for both negative and positive impact of radio-AGN outflows on star-fonnation. 
We therefore investigate the broadband spectral energy distributions (SEDs) of a sample of radio 
sources and try to decouple, by means of SED template fitting, the AGN and star-formation com- 
ponents. We use the cosmological parameters Hq = 71 kms^^Mpc^^, Q.m = 0.21 , and Q.\ = 0.73 
(from the first-year WMAP observations; Spergel et al. 2003). 

2. The AKARI North EcUptic Pole (NEP) field 

The AKARI space telescope (Murakami et al. 2007), using the InfraRed Camera (Onaka et al. 
2007) carried out observations of legacy fields, achieving greater sensitivity and resolution than 
previous instruments. One of these is the North Ecliptic Pole (NEP) field, which totals an area of 
~ 5.4 deg^ (Kim et al. 2012). With 9 spectral bands, ranging from 2.4jLtm to 24jLtm, the IRC con- 
tinuously covers the whole near- to mid-IR wavelength range, including the prominent wavelength 
gap (9-20/xm) that characterized Spitzer observations. 

For the NEP-Wide field exist also: deep GALEX observations (PI: Malkan), deep optical observa- 
tions (CFHT; Hwang et al. 2007, Maidanak; Jeon et al. 2010), near-IR observations (FLAMINGOS; 
Jeon et al., in prep.), and radio observations at 1.4 GHz (WSRT; White et al. 2010). The all-sky 
sur^vey of the WISE telescope (Wright et al. 2010) can provide us with additional data in the near- 
to mid-IR regime. Finally, several spectroscopic campaigns have also taken place (with WYIN, 
MMT, and Keck telescopes). In total more than 2000 spectroscopic redshifts are available in the 
NEP-Wide field (Shim et al., in prep.). 

3. Sample and method 

3.1 Radio-IR source cross-matching 

We cross-match the AKARI catalog of Kim et al. (2012) with the 1.5 GHz catalog of White 
et al. (2010). The original IR band-merged AKARI NEP-Wide catalog contains 114794 sources 
at an N2 {lAjXm) band AB magnitude limit of ~ 21 and a resolution of ~ 4 arcseconds. For the 
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cross-matching we use a sub-sample of 107504 near-IR detected sources (5a detection in N2 or 
N3 bands). The radio catalog covers an area of 1.7 deg^ with a beam size of 17 arcseconds and a 
sensitivity of 21 /i Jy/beam. The final catalog contains 462 radio sources at a 5a detection limit. 
To cross-match the two catalogs, we use the Poisson-probability based method of Downes et al. 
(1986) that takes into account both the proximity and apparent magnitude of an infrared neighbor 
source in order to calculate the probability of a true match. The final IR-radio catalog contains 425 
cross-matched sources. 

3.2 Photometric redshifts and radio-loudness 

We calculate photometric redshifts for the IR-radio cross- matched sources using the publicly 
available LePhare code (Arnouts et al. 1999, Ilbert et al. 2006). We use the neai-UV GALEX band 
as well as the full optical bands and near-IR bands, extending out to the W2 WISE band (4.6 /im). 
For the photometric redshift estimation we use the set of CFHT galactic SED templates from Ilbert 
et al. (2006), the PoUetta et al. (2007) AGN templates, and several stellar- template libraries. After 
excluding stars and sources for which LePhare failed to determine a photometric redshift, we are 
left with 266 sources with acceptable photometric redshifts. Given the limitations of photometric 
redshifts for the following we consider sources at z ^ 2.0. 

Originally radio-loudness, is defined as the ratio of the luminosities at 5 GHz and 4000A. 
Under that definition, radio-loud sources are those with Ri > 10. Given the available data for the 
NEP-Wide field, here we use an alternative definition of radio-loudness from Ivezic et al. (2002), 
utilizing the luminosities at 1.5 GHz and the CFHT i-band (7000A). Under that definition, we clas- 
sify sources with Rj > 2 as radio-loud, while for /?, < 1 a source is classified as radio-quiet. For the 
calculation of the radio-loudness we are using i-band fluxes corrected for Galactic extinction (from 
Schlegel et al. 1998). In Fig. |l| (left) we show the distribution of optical i-band and radio 1.5 GHz 
fluxes, together with constant radio-loudness lines. According to our definitions above, we have in 
total 70 radio-loud and 105 radio-quiet sources. 

3.3 SED fitting 

Using the full broadband SEDs (ranging from UV to mid-IR) we employ the method and tem- 
plates of Trichas et al. (2012) to fit sources simultaneously and additively with both an AGN and 
star-forming template. In this way, for each source the total bolometric luminosity, Lhoi, and the 
fractional bolometric contribution a of the AGN component are derived. An example of such a fit 
is shown in Fig. |I| (right). Through and eye evaluation of the fits we define a sub-sample of 156 
for which we have satisfactory a SED fit, of which 145 have redshifts z ^ 2.0. We use these for the 
following. 

4. Results 

4.1 AGN contribition vs. Li ^ghz 

In using the SED fitting results we can first of all see how the AGN and star-formation bolo- 
metric luminosity fractional contribution behaves as a function of both radio luminosity and radio- 
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Figure 1: Left: Extinction-corrected optical i-band flux and radio 1.5 GHz flux for the sample of 266 
IR-radio sources with acceptable photometric redshifts. The diagonal lines denote constant radio-loudness 
values, with the solid lines marking the limits for radio-loud and radio-quiet sources. Right: Example of a 
source SED and the total fit (blue) with both an AGN (green) and star-formation (red) contribution. 



loudness. The results are shown in Fig. § 
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Figure 2: AGN bolometric contribution as a function of 1 .5 GHz luminosity Li .sghz (left) and radio loudness 
Ri (right). The red solid line shows average values of AGN fractional bolometric luminosity contribution 
over Li sGHz bins. In the left panel with the vertical black dashed line is a transition luminosity taken from 
Mauch & Sadler (2007) (see text). In the right panel we also note the definitions of radio-loud and radio-quiet 
sources with vertical solid black lines. Horizontal error bars denote the respective bin width. 



Both for the highest radio luminosities and for the highest radio loudness values, sources have 
their bolometric luminosities dominated by an AGN component. These are the radio-AGN systems 
the star-formation properties of which we are interested in. The change between star-formation 
and AGN dominance^ in temis of bolometric luminosity agrees well with previous studies and lies 



We note that if a is the AGN fractional bolometric contribution, then the star-formation fractional bolometric 
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roughly ai-ound 10"*" erg/s or 10^-* W/Hz (e.g., Mauch & Sadler 2007). 

Several other points of interest arise from Fig. |[ Even at the lowest radio luminosities, where 
star-formation should dominate, we find on average ~ 40 ± 3% contribution to the energy output 
from an AGN component. Similar behavior is seen for the most radio-quiet sources. Conversely, 
for the highest radio-luminosities ~ 30 it 3% of the bolometric luminosity is contributed by star- 
formation. For the most radio-loud systems the star-formation contribution is further suppressed, 
with below ~ 20% fractional bolometric contribution. It should be noted that for the last bins of 
both radio luminosity and radio-loudness, we observed a change in the aforementioned trends. This 
is most probably an effect of low-number statistics (2 and 7 sources respectively contained in each 
of these two bins) combined with possible photometric redshift effects. 




Figure 3: Ratio between observed and expected star-formation rates (see text) as a function of Li ^ghz (left) 
and as a function of radio-loudness Rj (right). We plot average values over luminosity and radio-loudness 
bins, respectively, for AGN dominated sources. Numbers denote the AGN dominated average redshifts and 
number of sources for each respective bin. The horizontal line separates the negative (below) and positive 
(above) feedback space. The vertical lines are the same as in Fig. 2. 



5. Star-formation quenching? 

We use the star-fomiation component far-IR luminosity ^s-iooo/im from the SED fitting to cal- 
culate the SFR for the IR-radio sample using the relation from Kennicutt (1998). Although this is 
a theoretically derived relation, empirical calibrations of the relation between SFRs and Lg-iooo^m 
exist in the literature. These are usually within 30% of the relation from Kennicutt (1998) and 
therefore broadly agree with each other. 

In Fig. ^ the ratio between observed and expected SFRs per radio luminosity (left) and radio- 
loudness bins (right) is presented for the AGN dominated objects^. By using the empirical relation 
from Elbaz et al. (201 1) that connects specific star-formation rate (sSFR) with z, we can calculate 

contribution is, by definition, 1 — a. 

^We define AGN dominated objects such tliat their AGN fractional bolometric contribution is larger than 55%. 
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the expected SFRs for a passively evolving galaxy at the given redshift of each bin and assuming a 
stellar mass of = 10^^ '^Msoi (mean stellai^ mass of radio-loud AGN, e.g., Best et al. 2005). We 
can thus calculate the ratios plotted in Fig. |[ 

When looking at objects with luminosities around ~ 10^^ erg/s, we see that with an average redshift 
of ~ 1, they exhibit SFR below the expected levels. The same behavior is also seen when consid- 
ering radio-loudness, with sources at 7?, '--^ 2.5 exhibiting lower SFRs than expected given their 
average redshift. However, if we look at the last bin of both radio luminosity and radio-loudness, 
objects appear to show the expected amount of star-formation. 



6. Conclusions 



In conclusion, we see hints for lower SFRs than expected for radio-AGN dominated systems. 
However, if we consider the statistical errors together with the intrinsic scatter of the assumed sSFR 
to z empirical relation, we do not see a strong signature of negative feedback in radio-loud AGN. 
Conversely, our data clearly demonstrate the absence of positive feedback in radio-loud AGN. Es- 
timation of individual stellar masses for each of our objects will allow us to differentiate between 
a negative or a non-existent feedback signature in these sources. 
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